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Already more that 40 years ago, it has been suggested that because of the enormous mass densities 
in the cores of neutron stars, the hadrons in the centers of neutron stars may undergo a phase 
transition to deconfined quark matter. In this picture, neutron stars could contain cores made 
of pure (up, down, strange) quark matter which are surrounded by a mixed phase of quarks and 
hadrons. More than that, because of the competition between the Coulomb and the surface energies 
associated with the positively charged regions of nuclear matter and negatively charged regions of 
quark matter, the mixed phase may develop geometrical structures, similarly to what is expected 
of the sub-nuclear liquid-gas phase transition. In this paper we restrict ourselves to considering the 
formation of rare phase blobs in the mixed quark-hadron phase. The influence of rare phase blobs 
on the thermal and transport properties of neutron star matter is investigated. The total specific 
heat, cv, thermal conductivity, k, and electron-blob Bremsstrahlung neutrino emissivities, £„,br, of 
quark-hybrid matter are computed and the results are compared with the associated thermal and 
transport properties of standard neutron star matter. Our results show that the contribution of rare 
phase blobs to the specific heat is negligibly small. This is different for the neutrino emissivity from 
electron-blob Bremsstrahlung scattering, which turns out to be of the same order of magnitude as 
the total contributions from other Bremsstrahlung processes for temperatures below about 10 8 K. 

PACS numbers: 21.65.Qr; 26.60.Gj; 97.10.Cv; 97.60.Jd 



I. INTRODUCTION 

Already many decades ago, it has been suggested that, 
because of the extreme densities reached in the cores of 
neutron stars, neutrons and protons may transform to 
quark matter in the cores of such objects [H-Q- Quark 
matter could thus exist as a permanent component of 
matter in the ultra-dense centers of neutron stars 

. If the dense interior of a neutron star is indeed con- 
verted to quark matter, it must be three-flavor quark 
matter since it has lower energy than two-flavor quark 
matter. And just as for the hyperon content of neu- 
tron stars, strangeness is not conserved on macroscopic 
time scales, which allows neutron stars to convert con- 
fined hadronic matter to three-flavor quark matter until 
equilibrium brings this process to a halt. 

As first realized by Glendenning the presence of 
quark matter enables the hadronic regions of the mixed 
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phase to arrange to be more isospin symmetric than in 
the pure phase by transferring charge to the quark phase 
in equilibrium with it. The symmetry energy will be 
lowered thereby at only a small cost in rearranging the 
quark Fermi surfaces. The electrons play only a minor 
role when neutrality can be achieved among the baryon- 
charge carrying particles. The stellar implication of this 
charge rearrangement is that the mixed phase region of 
the star will have positively charged regions of nuclear 
matter and negatively charged regions of quark matter. 

Because of the competition between the Coulomb 
and the surface energies associated with the positively 
charged regions of nuclear matter and negatively charged 
regions of quark matter, the mixed phase may develop 
geometrical structures (see Fig. [I}, similarly as it is ex- 
pected of the subnuclear liquid-gas phase transition [TTr - 
[l3j . This competition establishes the shapes, sizes, and 
spacings of the rare phase in the background of the other 
in order to minimize the lattice energy [3, H3, E3] • 

The change in energy accompanied by developing such 
geometrical structures is likely to be very small in com- 
parison with the volume energy [Io|, and, thus, 
may not much affect the global properties of a neutron 
star. However, the geometrical structure of the mixed 
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TABLE I. Masses tub, electric charges Qb, spin Jb, and third 
component of isospin Ig of the baryons B included in the 
lagrangian of Eq. fTJ) [3, [{| . 
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FIG. 1. Schematic illustration of possible geometrical struc- 
tures in the quark-hadron mixed phase of neutron stars. The 
structures may form because of the competition between 
the Coulomb and the surface energies associated with the 
positively charged regions of nuclear matter and negatively 
charged regions of quark matter. 



phase may be very important for irregularities (glitches) 
in the timing structure of pulsar spin-down as well as 
for the thermal and transport properties of neutron stars 

BE El- 

To calculate the neutrino-pair bremsstrahlung rates 
and thermal properties, we follow the method described 
in [TH and 19], which is commonly used for the calcula- 
tion of the neutrino emissivity and thermal conductivity 
in the crusts of neutron stars. These authors considered 
contributions from electron-phonon scattering and Bragg 
diffraction (the static-lattice contribution). Furthermore, 
multi-phonon processes and electron band structure ef- 
fects are incorporated to obtain more realistic scattering 
rates and a better connection between the solid and the 
liquid gas phase. Instead of adopting the analytic fits 
provided in [lj| and [ijj], here we re-calculate the scat- 
tering rates from phonon sums using the method of pfjl ] . 
There are two main reasons for this. The first being 
that, for the crust, the total ion charge is balanced by 
the total electron charge. This will be different for the 
mixed quark-hadron phase in the core of a neutron star, 
since electric charge neutrality is established between the 
electric charges of the rare phase, the dominant phase, 
and the leptons which are present in both the rare and 
the dominant phase. The simple relation n e = Zni be- 
tween electron density and ion density, used to derive the 
crustal fit formula in [18L , can therefore not be used to 
study the quark-hadron Coulomb lattice structure in the 
core a neutron star. The second reason concerns the elec- 
tric charge numbers themselves. For mixed phase blobs, 
they can easily exceed Z ~ 10 3 , as will be shown in § 111 CI 
Charge numbers that high are obviously not reached in 
the crustal regimes of neutron stars where there is 
usually no need to consider atomic nuclei with charges 
much larger than Z > 56. 

The paper is organized as follows. In Section [Til we 



briefly discuss the modeling of the mixed quark-hadron 
phase in the cores of neutron stars and the equations 
of state of confined hadronic and quark matter used in 
this work. In Section IIII1 we summarize the formalism 
for calculating the neutrino-pair Bremsstrahlung emis- 
sivity and the thermal conductivity of rare phase blobs 
immersed in hadronic matter. The results are presented 
in Section Hvl 



II. MODELING OF THE MIXED 
QUARK-HADRON PHASE IN NEUTRON STARS 

A. Hadronic matter 

To compute the particle compositions of the cores of 
standard neutron stars, that is, neutron stars without 
deconfined quark degrees of freedom, we choose a rela- 
tivistic lagrangian of the following type @, Q , 

B 

- (m N - ffo-cr)]^ + i(9 M 0-d'V - to 2 <j 2 ) 

- ^b a m N (g a a) 3 - ^c a {g a a) A - \j^ v ^ v 

+ + \m% ■ r (1) 

X—e~ 

where the sum over B sums the baryon species listed in 
Table HI The sum over A accounts for the presence of 
relativistic electrons and muons in neutron star matter. 
Their masses are m e = 0.511 MeV and = 105 MeV. 
The quantities g p , g a , and g^ are meson-baryon coupling 
constants of er, l</', and p 1 " mesons. Non-linear a- meson 
self-interactions are taken into account in Eq. ([IJ via the 
terms proportional to b a and c a . The quantities r and 7 M 
denote isospin vectors and Dirac matrices, respectively, 
and f9 M = d/dx^ @, |||. We have solved the equations 
of motion for the baryon and meson fields, which follow 
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TABLE II. Relativistic mean-field parametrizations used in 
this work. 





Parametrizations 


Coupling constants 


HV 


G300 




8.7982 


9.1373 




9.1826 


8.6324 


9p 


9.7145 


8.3029 


bo 


0.00414 


0.003305 




0.00716 


0.01529 



from Eq. ([T]), in the framework of the relativistic mean- 
field approximation [3, E|, where the fields a, ui, p are 
approximated by their respective mean-field expectation 
values a = (a), u> = (u), and p = (po3)- Two popular 
jarametrizations labeled G300 and HV have been used 
3, IS [HUH ■ Tneir parameters are summarized in Table 

Neutron star matter is characterized by the conserva- 
tion of two charges, electric and baryonic. This feature 
leads to the chemical equilibrium condition 



Pi Bip n QiPe j 



(2) 



with p n and p e are the chemical potentials of neutrons 
and electrons. The quantities Bi and Qi stand for 
the baryon number and the electric charge of particles 
(mesons and baryons) of type i. Equation (J5|) greatly 
simplifies the mathematical analysis, since only knowl- 
edge of two independent chemical potentials, p n and p e , 
is necessary. The latter are given by 



Pb = 9uU + 9 P Pa?,I'B + 
px = 



K B 



l B ' 



•x • 



(3) 



where m* B = raj - g a B& denote the effective medium- 
modified baryon masses, ks and k\ are the Fermi mo- 
menta of baryons and leptons, respectively, and Ig is 
the third component of the isospin vector of a baryon 
of type B. Finally, aside from chemical equilibrium, the 
condition of electric charge neutrality is also of critical 
importance for the composition of neutron star matter. 
It is given by 



Y,Qi(2j B + r 



k 3 

K B 

6tt 2 



3tt 2 



0. 



(4) 



Figure [2] shows the baryon-lcpton compositions of neu- 
tron star matter computed from Eq. ([T]) for the relativis- 
tic mean-field approximation. The quantity pi in Fig. [5] 
stands for the individual number densities of baryons, 



p B = (2J B + 1) k 3 B /3iT 3 . 
whose total number density is given by 

P b = ^2 P B ■ 

B 



(5) 



(6) 
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FIG. 2. (Color online) Sample baryon-lepton compositions, 
pi/pb, of neutron star matter computed for G300 (upper 
panel) and HV (lower panel). 



The individual number densities of electrons and muons 
(i = e~ , p~) are given by 



P i = 2k 3 /3n 3 



(7) 



The total energy density and pressure of the matter, 
shown in Fig. [21 follow from 



+ \( m ^) 2 + -j{m p pf 

(■ks 



and 



PH = -^bm N {g (J a) 3 - ^c{g a a) A - i(m CT cr) 2 



(8) 



k 2 dk\ k 2 + m 2 



(9) 
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B. Quark matter 

To model quark matter, we use the MIT bag model. 
Up (u) and down (d) quarks are treated as massless parti- 
cles while the strange quark (s) mass is assigned a valued 
of m s = 200 MeV. First-order perturbative corrections in 
the strong interaction coupling constant a are taken into 
account |24T - l27j . The Landau potentials of up and down 
quarks are then given by 



47T 2 \ 7T 



4 



4;r 2 \ ' 7T J 



(10) 



(11) 



while for strange quark we have 



n„ = — 



2 s 



(m*V 



Ms ~ m l ~ mj(n s ,m a ) 



2a 

7T 

2(^ -m^ + Sm* In 2 ^ 

6 In — ( fi s m 2 8 ^ p? s -m? s - mif(fi„, m s ) ) } 



(12) 



where f(p,m) = ln((y« 4- y y« 2 — m 2 )/m), and ct is a 
renormalization constant whose value is of the order of 
the chemical potentials [261 ]. In this article we take 
ct = 300 MeV. The Landau potentials of electrons and 
muons are given by 




(13) 



— m,, In 

1 M I 



- J- m 2 K 



(14) 



The condition of chemical equilibrium leads to 

Pd = Ms = Mu + Me = M« + Mm • (15) 

The partial baryon number densities of the particles is 
obtained from (i — u, d, s, e — , fi~) 

Pi = -dQ,i/ d^i, (16) 

and the total energy density and pressure of quark matter 
follows from 



and 



(17) 



(18) 



where i? denotes the bag constant. For a,m s — > one 
recovers from Eqs. (IT7|) and (fT5|) the standard equation of 
state of a massless relativistic quark gas, P = (e — 4B)/3. 



C. Geometric structures in the mixed 
quark-hadron phase 

To determine the possible geometric structures in the 
mixed phase of quarks and hadrons, we use the Gibbs 
condition 



Pif(//n,Me,{0}) = Pq(Pu, Me 



(19) 



for phase equilibrium between hadronic matter and quark 
matter The quantity {</)} in Eq. (Tl9l) stands collec- 
tively for the field variables (ct, Q, p) and Fermi mo- 
menta (fes, &a) that characterize a solution to the equa- 
tions of confined hadronic matter. We use the symbol 
X = Vq/V to denote the volume proportion of quark 
matter, Vq, in the unknown volume V. By definition, \ 
then varies between and 1 , depending on how much con- 
fined hadronic matter has been converted to quark mat- 
ter. Equation (fT9]) is to be supplemented with the con- 
dition of baryon charge conservation and electric charge 
conservation. The global conservation of baryon charge 
is expressed as 

Pb = XPQ{Pn,Pe) + (1 - X) PH{Pn,Pe, {<t>}) , (20) 

and the global neutrality of electric charge is given by 

= X <7Q(Mn,Me) + (1 - X) qH(Pn, Me, {</>}) • (21) 

In Figs. [3] and [4] we show sample compositions of neu- 
tron star matter computed for four different parameter 
sets (see Table IIIII) which allow for the presence of a 
mixed phase of quark and hadrons. As can be seen from 
these figures, the mixed quark-hadron phases obtained 
for the G300I and G300II parametrizations exists over a 
broader density range than it it the case for HV1 and 
HV2. This feature has it origin in the density depen- 
dence of p n , which increases faster with density for the 
HV1 and HV2 models, rendering these equations of state 
stiffcr than the G300I and G300II models for the equa- 
tion of state. 

A larger value for the bag constant, B, reduces the 
pressure of the quark phase, which too leads to a broader 
mixed phase region. It also follows from Figs. |3] and 0] 
that the number of electrons as well as the number of 
muons drops quickly with density, since charge neutrality 
is achieved chiefly among the quarks themselves. For 
values of the strong interaction coupling constant a > 
0.25 electrons disappear from the matter and positrons 
tend to emerge in the mixed phase. We therefore consider 
only a < 0.2. 

In Figure[S]we show the masses of neutron stars, whose 
compositions are given in Figs. [3] and [4] The underlying 
equations of state are given in Eqs. ©, ©, (jTT)) and 
(|T£|). and the parameter sets are listed in Table [TTT1 The 
Baym-Pethick-Sutherland (BPS) model for the equation 
of state has been used to model the crusts of these neu- 
tron stars [28]. The maximum masses of the quark- 
hybrid stars computed for HV1 and HV2 are 1.47 M Q 
and 1.61 M Q , respectively. For G300I and G300II we 
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FIG. 3. (Color online) Particle composition, Pi/pt, of neutron 
star matter with mixed quark-hadron phase. The parameter 
sets are HV1 and HV2 (see Table HIT|| . 



TABLE III. Parameter sets of the quark-hybrid star models 
of this work. B denotes the bag constant, a is the strong 
interaction coupling constant. 



Label 


Hadronic Phase 




Quark Phase 






HV1 


HV 


B = 


110 MeV fm" 3 , 


a 


= 0.2 


HV2 


HV 


B = 


160 MeV fm~ 3 , 


a 


= 0.1 


G300I 


G300 


B = 


110 MeV fin" 3 , 


a 


= 0.2 


G300II 


G300 


B = 


160 MeV for 3 , 


a 


= 0.1 



obtain maximum masses of 1.58 M Q and 1.69 M Q , re- 
spectively. These values are too low to accommodate the 
recently discovered heavy pulsar PSR J1614-2230, whose 
mass is M = 1.97 ± 0.04 M [29|. One possible explana- 
tion could be that the high rotation rate of this neutron 
star prevents the hadrons in the core of this neutron star 
from transforming to quark-hadron matter. This neu- 
tron star would thus be made of confined hadronic mat- 
ter, whose equation of state is stiffer than the equation 
of state of quark-hybrid matter, supporting stars close to 
2 M.0 for the two models considered in this paper (curves 
labeled HV and G300 in Fig. [5]). As found in (30|, massive 
(~ 2 M Q ) non- rotating neutron stars with extended re- 
gions of deconfined quarks and hadrons are comfortably 
obtained in the framework of the nonlocal SU(3) Nambu- 
Jona Lasinio model. This model is not considered in this 
paper, however, since our results are largely independent 
of the particular microscopic many-body model chosen 
to determine the equation of state of bulk quark-hadron 
matter. 
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FIG. 4. (Color online) Same as Fig. O but for parameter sets 
G300I and G300II (see Table |HT)|. 

As is pointed out in Q the isospin restoring force can 
exploit degrees of freedom made available by relaxing the 
requirement of strict local charge neutrality to neutrality 
on larger scales and form a positively charged hadronic 
matter region with lower isospin asymmetry energy and 
a negatively charged quark matter region. The compe- 
tition between the Coulomb interaction and the surface 
energy will result in a crystalline lattice of the rare phase 
in the dominant phase. The size and spacing of the 
crystalline lattice is determined by minimizing the to- 
tal energy. The situation is similar to the atomic nuclei 
immersed in a relativistic electron gas in the crusts of 
neutron stars. 

Depending on density, the embedding of the rare phase 
in the dominant phase can lead to different geomet- 
ric structures, including spherical blobs, rods and slabs. 
In what follows, we restrict ourselves to the discussion 
of spherical blobs. For an electrically charge neutral 
Wigner-Seitz cell, with spherical blob of radius r^, the 
Coulomb and surface energy density can be expressed as 

e c = 2na e [q H { X ) ~ q Q (x)} 2 r\ X h{x) , (22) 
es = 3 X a(x)/r & , (23) 

where a e = 1/137 is the fine structure constant, x = 
min(x, 1— x) is the volume fraction of rare phase, fs(x) = 
(x— 3x 1 / 3 +2)/5 is the function fdix) for d = 3 which arise 
from calculating the electrostatic binding energy of the 
cell @. Due to theoretical difficulties it is very hard to 
estimate the surface tension a, here we follow p} and take 
a gross approximation expression for the surface tension 
first proposed by Gibbs [31| , where the surface energy is 
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FIG. 5. (Color online) Mass-central density (top panel) and 
mass-radius (bottom panel) relationships of neutron stars for 
the parameter sets listed in Table Hill The solid dots denote 
the most massive neutron star model of each sequence. 



FIG. 6. (Color online) Radii of spherical blobs, r, and Wigner- 
Seitz cells, a, as a function of quark volume fraction, Xi f° r 
different surface tensions, r\. The parameter sets are listed in 
Table [ml 



proportional to the difference of the energy densities of 
the two phases, 

a(x)=vL[e Q (x)-e H ( X )} , (24) 

where r\ should be on the order of 77 ~ 0(1), and we 
take L = 1 fm. Three different values for the constant 
77, 0.5, 1, and 2, are used in our calculations in order to 
investigate effects caused by uncertainties in the value of 
the surface tension. 

Since ec oc r\ and es oc r^ 1 , it is possible to minimiz- 
ing the total energy ec + es at fixed \ which leads to 
an equilibrium radius of the rare phase of blobs inside of 
Wigner-Seitz cells, 



\^a e [q H {x) - mix)] h(x) J 

The radii of spherical blobs of the rare phase, rt,, and 
the radii of Wigner-Seitz cells, a, as a function of the 
quark volume fraction x are shown in Fig. [5] The radii of 
spherical blobs is in the range of 10 to 30 fm. To compare 
our situation to the crust we also calculate the charge Z 
and mass number A = mb/m u of the rare phase blobs, 
with mf, being the blob mass. In addition we define an 
effective electric charge number given by 

£eff = -, (26) 

where n& is number density of the spherical blobs. The 
values of Z , A and Z c g as a function of the quark volume 



fraction are shown in Figs. [8] and [9] As can be seen, the 
charge and mass numbers are typically one to two orders 
of magnitude greater than that of the heaviest stable nu- 
clei, but due to a dramatic drop in electron density the 
value of Z c ff can fall near and below Z e g ~ 10 as x — >■ 1. 
The discontinuities of the curves at x — 0-5 in Figure [7] 
is due to the differences in the hadronic and quark phase 
densities at x = 0.5. It is also worth noting that near 
the edges of the mixed phase region (x —> 0, 1 or x —> 0) 
the volume density of blobs rib = 3a;/ (47rrf ) vanishes, but 
since /a(0) = 2/5, the blob radius r& approaches a con- 
stant rb — > rb(x = 0). Therefore Z c g = n e /rib diverges at 
the edge of the mixed phase region. 



III. THERMAL AND TRANSPORT 
PROPERTIES OF QUARK-HADRON PHASE 

Next, we turn our interest to the calculation of the 
thermal and transport properties of a mixed phase of 
quarks and hadrons. Knowledge of these properties is of 
key importance in order to carry out thermal evolution 
simulations of neutron stars with hypothetical quark- 
hadron cores and to determine possible astrophysical sig- 
natures hinting at the existence of such matter inside of 
neutron stars. Our focus here is on exploring the impact 
of rare phase blobs on the following properties: specific 
heat cy, neutrino emissivity e„, and thermal conductivity 

K. 

From the outset, one might expect that, because of the 
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FIG. 7. (Color online) Mass number, A = mt/m u , of spheri- 
cal blobs of rare phase as a function of quark volume fraction, 
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FIG. 8. (Color online) Electric charge number, Z, of spherical 
blobs of rare phase as a function of quark volume fraction, \- 



geometric structures in the mixed phase, new degrees of 
freedom are introduced to the system, which may store 
additional thermal energy and hence would increase the 
specific heat. Due to the scattering of degenerate elec- 
trons and rare phase blobs in the mixed phase, there 
will also be an additional contribution to the neutrino 
Bremsstrahlung and an additional microscopic entropy 



FIG. 9. (Color online) Effective charge number, Z e s, of spher- 
ical blobs of rare phase as a function of quark volume fraction. 
X- The five horizontal lines at the bottom of each panel cor- 
respond to Z^ff ' min for n = 1 , 2, 3, 4, 5 (see Section [III Bp 



production rate which reduces the total thermal conduc- 
tivity. 



A. Basic physical quantities 

We hrst introduce some of the basic physical quanti- 
ties and functions used in the calculation of the thermal 
and transport properties of the quark-hadron phase. The 
state of the mixed phase with rare phase blobs is deter- 
mined by the ion-coupling parameter (l8| . 



r = 



Z 2 e 2 
aksT 



(27) 



where a — (3/(47TOb)) 1 / 3 is the Wigner-Seitz cell ra- 
dius, which is related to the spherical blob radius by 
x = (ri,/a) 3 . The quantity T rn = 172 corresponds to 
the melting point below which a classical one-component 
Coulomb crystal (T < r m ) becomes a Coulomb liquid 
(1 < r < T m ) [32]. Therefore the melting temperature 
is given by T m = Z 2 e 2 /(akB^m)- Another important 
quantity is the plasma temperature, 



Tn 



(28) 



where uj p — (4irZ 2 e 2 rib/mb) 1 / 2 denotes the plasma fre- 
quency of the spherical blobs, with mj the mass of the 
spherical blobs. For later, it is convenient to introduce 
the temperature in units of the plasma temperature, 
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t p = T/Tp. For temperatures t p < 1/8 the vibration of a 
Coulomb crystal must be treated quantum mechanically 
18]. 

Besides the Coulomb interaction, there are three effects 
which must be taken into account to describe the inter- 
action between electrons and rare phase blobs. These are 
the screening of electrons, the shape of the blob and the 
effect of thermal vibration. The Fourier transform of the 
effective electron-blob interaction is given by 18] 



= ^epzFjq) e _ w{q) 



(29) 



where pz is the blob charge per unit volume pz = X1Q 
for x < 1/2, and p z = (1 - x)qh for x > 1/2, The 
quantity e(q) in Eq. (|29l) is the static longitudinal dielec- 
tric factor adopted from t 33] and F(q) is the form factor 
of a blob. For simplicity, we assume a uniform distri- 
bution of the electric charge in the rare phase and use 
F(q) = (3/(qr) 3 )[sm(qr) - qr cos(qr)] QJ]. 

Thermal vibrations of rare phase blobs are taken into 
account via the Debye- Waller factor, 



W(q) = 



hq 2 I coth(hu s /(2k B T)) 
4mb \ w s 



(30) 



where uj s is the phonon frequency. Here (. . .) p h denotes 
the average over the phonon wave vectors and polar- 
izations. Throughout this article, we use the method 
of [20j to compute phonon sums. It is assumed that 
there are three polarizations of phonons: two trans- 
verse modes with linear dispersion relations = a,fc 
(i = 1,2), and one longitudinal mode which is deter- 
mined through Kohn's sum rule uj 2 1 + lu 2 2 + oj 2 = w 2 , 
where lj p denotes the plasma frequency. The two param- 
eters a± and a 2 are determined by fitting the frequency 
moments u n = (u/ l ) p h of the specified lattice type. For 
a bec lattice, the frequency moments u_i = 2.7990 and 
m_2 = 12.998 are well known (34j, and are used to ob- 
tain ai — 0.58273, a 2 — 0.32296. These parameters 
also produce first and fourth moments: p\ = 0.51106, 
P4 = 0.201946 which are quite close to exact the values 
pi = 0.51139 and p4 = 0.203076 34]. In the special case 
of the Debye- Waller factor, the phonon sum can be fitted 
very well by the following analytic formula p35j . 



W(q) 



= «o / _q_\ f i 



2 \2k F J \2 



(31) 



where fcp is electron Fermi wave number and ao is a 
constant given in [l8|. The latter can be rewritten as 



ao _ 3 2 / 3 7r 5 / 6 ^/ 2 nT 
2 Oie^c 1 / 2 Zrn^ 2 n^ 



(32) 



With the help of Eq. (|3"2")l the Debye- Waller factor can be 
written as W(y 2 a,t p ), where y — q/(2k F ). 



B. Neutrino Bremsstrahlung emissivity 

Next we turn to the calculation of the scattering 
of electron off rare phase blobs (that is, electron- 
blob Bremsstrahlung), which leads to the generation of 
neutrino-anti-neutrino pairs according to the reaction 
e~ + (Z, A) — > e~ + (Z, A)+i>+D. The associated neutrino 
Bremsstrahlung emissivity can be written as t 36] 



tTrG 2 Z 2 e 4 C 2 



■{k B Tfn b L.. 



(33) 



where rib is the number density of the rare phase blobs 
and G F = 1.436 x 10~ 49 er g cm 3 is the Fermi weak cou- 
pling constant, C\ « 1.675 [To ], and L is a dimensionless 
function given by 



L = L 



^liq 



(34) 



where L p h accounts for the scattering of electrons off the 
phonons of the Coulomb crystal of rare phase blobs, L s \ 
accounts for Bragg scattering between electrons and the 
static Coulomb crystal lattice, and Ln q is for liquid phase. 
In the liquid phase the general expression is obtained 
through a variational approach in the Born approxima- 
tion l36Tj38ll, 



2y 2 



1-y 2 



Iny), (35) 



where y = q/(2kp)- We follow [18] for the choice of the 
ion-ion structure factor S{q), which was fitted in (39l.l40|. 
For a solid phase the phonon contribution is mainly given 
by umklapp processes which, in Born approximation, can 
be written as fTHH 



i P h = dy 



S eB (y 2 a,t p )\F(y)\ 



Vo 



yHvW 



2y 2 



1 



y 



■lny 



(36) 

where the lower integration limit yo = (4Z ff) 1 / 3 ex- 
cludes the low-momentum transfers in which the umk- 
lapp processes are forbidden 18], and the effective static 
structure factor S c g is obtained from the summation of 
multiphonon diagrams |18l |4l| . For the parameter sets 
chosen in this calculation we always have Z g > 1/4 so 
that yo < 1 . 5 c ff can be written into the form of a rapidly 
decreasing integral |18| . 

S cS (y 2 a,t p ) = 189 ( -X e- 2W ^ & ^ 



< / d£ 



1 - 40<f + 80e 



(l+4e 2 ) 5 cosh 2 «) 



where 



(t) 



) sinh ( ; 



(37) 



(38) 



P h 
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with oj = LOs/cdp denoting the phonon frequency in units 
of the plasma frequency. Similar to the analysis in [l8[ , 
there exist approximate expressions for .S^ff for the lim- 
iting cases where t p <^ 1 and t p 3> 1, which is discussed 
next. For this purpose, we write the cc-independent part 
of $(£,x,tp)/t p as 



1 



p \uj sinh 



(2*,,) / P h 



(39) 



The limits t p <C 1 and t p 3> 1, ip{^,tp) can then be 
replaced by 



2 ( 4 °) 



where 



^o,t p ) - F(t p )e 1 tp»i, 



t P -!-o i/}(0,t p ) 



(41) 



is computed numerically. It is a rapidly decaying function 
of £ and is negligibly small for £ > 2. The function 
F(t p ) = (u)/[2t p sinh(uj/2t p )})ph is computed numerically 
for 1 < t p < 10 2 . Asymptotically, F(t p ) has the form 



(42) 




*p ^ 10 2 ■ 



wi p sinh I 



21, 



(43) 



P h 



can be calculated numerically for a broad range of t p 
values, and can be shown to behave as 



^(o,t P ) 



2u 



t P >l, 



2 12tf ' "P 

^ 2 (^ + ^K, t p <io-\ 



(44) 



With the aid of Eqs. (|3"9")l to (|44p . we can now derive 
low and high temperature limits of the effective structure 
factor S e g of Eq. ([37| . For t p <C 1 one obtains 



S cff (x, *„) = 189 0) e- 2M ^<^G o ffK^(0, t p )) , 

(45) 



while for t p ^ 1 



-2W(a:,t p )+a;tpi/>(a:,tp) 



SeffOMp) = 189 I - I < 

x H eS (xt p F(t p )) - e - 2W ^ . (46) 
Here we have defined 

1 - 40£ 2 + 80£ 4 



G c ff(a) 



(1 +4£ 2 ) 5 cosh 2 (7r£) 



x (e°*«> - l) , 



(47) 



TABLE IV. Lower bound values Z^^. 



7(1) 

^cff.min 


7 (2) 
^cff.min 


7(3) 
^cff.min 


7 (4) 
eff ,min 


7 (5) 


V2tt/3 


4tt/3 




8\/27r/3 


5^71/3 



which obeys 



Geff(a) 



181440 



a > 10 2 : 
a, a<0.1. 



(48) 



For 0.1 < a < 10 2 the value of G c ff is obtained numeri- 
cally. The quantity H e g in Eq. ([4T>)) is defined as 

/■°° 1 - 40£ 2 + 80£ 4 „ P % , , 

iJeff(a) = / 9 e" a? . (49 

V 7 Jo (l + 4£ 2 )5cosh 2 K) V 7 

Asymptotically, 7? e ff behaves as 

Hfr(a) ~ fi|9(f) 5 + 0.003690a, a < 0.1 , 

^ ff W-| (1/2)7r l/2 a -l/2 5 a > 1Q 3 5 i 5 °) 

but its values for 0.1 < a < 10 3 need to be computed 
numerically. 

To determined the contributions of the static 
lattice contribution (Bragg diffraction) to neutrino 
Bremsstrahlung, we follow [18], who considered band 
structure effects. We begin with defining the dimension- 
less factor L s \, 



\2Z c fr — ' 

K^0 



\<K)\ 



(51) 

where the sum is over K values below the electron 
Fermi surface, i.e. yn — K/(2kp) < 1 and ty = 
\V K \y/l-y 2 /(k B T) where V K is given in Eq. (gU). The 
integral I(y,ty) in Eq. (IBTj) has been fitted analytically 
over a wide range of ty and y values in (l8| . 

In our case, which deals with a solid phase of rare phase 
blobs immersed in hadronic matter, the summation of 
L s \ sometimes consists only of a few terms. In the case 
of a bec lattice, the condition for the summation to have 
at least n terms (not counting multiplicity) is yK n < 1, 
which translates to a lower limit for Z e rr > Z„. • . The 

— err, mm 

values of zj£' ' • for n — 1 through n = 5 are shown in 



Table EVJ 

Therefore, since Z c g = Z in a lattice of nuclei im- 
mersed in an electron gas, L s \ is guaranteed to have 
a handful of terms. However, as we have seen in Sec- 
tion III CI the electron density will drop drastically with 
increasing baryon number density in the mixed phase re- 
gion if x ^ 0.5 so that Z e g drops correspondingly (see 
Fig. [5] where the first few Z^ min are shown). 



As will be shown in Section [TV] for low temperatures 
of T < 10 8 K and a quark volume fraction \ <L 0.5 the 
summation consist only of a few terms, and the contri- 
bution of the static lattice oscillates vividly as a function 
of x (see Figs. E] and [H) 
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C. Thermal conductivity 

To calculate the thermal conductivity of quark-hybrid 
matter, we closely follow the formalism outlined in (l9j . 
In general, the thermal conductivity of degenerate elec- 
trons is expressed in terms of an effective collision fre- 
quency v K , 42] 



irkgTn e 



(52) 



which, in turn, can be expressed in terms of dimensionless 
Coulomb logarithms A K [ilj . 



4irZ 2 e i n b 

PpVF 



(53) 



where n& is the number density of rare phase blobs in the 
mixed phase, and pf and Vf are the Fermi momentum 
and Fermi velocity of rclativistic electrons. For a solid 
phase, the Coulomb logarithms are calculated variation- 
ally in Born approximation [l9l [42j , 



A 



K. solid 



dyS K (y)- 

yo 



y\<y)\ 2 



i - 



i 



(54) 



Here, x r = pF/{jn e c) is a relativistic parameter and yo, 
F, and e have the same meaning as in Section UlI Al The 
quantity S K stands for the effective static structure factor 
for thermal conductivity, [ljl |4l| 



S K (y) = S a (y)+[^--)SS K ( y ), 



(55) 



where y — q/{2kp) and S a {y) are effective structure fac- 
tors used to calculate the electric conductivity a [1 91 ] - 

Instead of using the asymptotic expressions and fitted 
formulas for A K provided in [19j . which cover 10" 3 < t p < 
10 and < ay 2 < 0.15 (a is given by Eq. ([32]) ). here we 
fully calculate their values and derive their asymptotic 
behaviors. This covers a wider range of t p and ay 2 val- 
ues. For thispurpose, we rewrite the relevant integrals 
given in (l9l l4lj in a form which is similar to S e g in Sec- 
tion HUH i.e., 



(y 2 a, t p ) = / 5— (e^*^ l) 

Jo cosh irt; v ' 



x e 



-2W(y 2 a,t p ) 



SS K (y 2 a,t p ) = d£ 



.2tt(1 - 2sinh 2 tt£) 

COSh 4 71"£ 



x e 



(56) 



where i/)(£,t p ) denotes the phonon sum function already 
defined in Eq. (f39|) . Since Bragg diffraction does not 
contribute to the thermal conductivity, there is no coun- 
terpart to L s \, 



Similar to Section MI B[ asymptotic expressions for 
tj){g,tp) may be used to derive the high and low tem- 
perature limits for S a and 8S K . One then obtains 

S a (x,t p ) ^> G a {xt p ^{Q,t p ))e- 2Wi - x ^ 

(e xt ^ Q ^H a (xt p F(t p )) - l) e~ 2W ^^ , 

SS K (x,t p ) Gs4^(0,t P ))e- 2W ^ 

> H SK {xt p F{t p ))e xt ^°' t ^-' 2W{ - :c ' t ^ , (57) 



where -F(ip) is given in Section IIII Bl and the functions 
Ga.Sn an d H<t.&k a re defined as 



G a (a) = 


/ 

Jo 


Gsn(a) = 


p oo 

/ d£ 

Jo 


H a {a) = 


/ d£ 

Jo 


H 5K (a) = 


/>oo 

/ dH 

Jo 



cosh" 7r£ V 
2tt(1 - 2sinh 2 7r£) a 

COSh 7r£ 



V>(£) 



cosh 7r£ 



2vr(l - 2sinh 2 7r£) _ Q? 2 
cosh 4 7r£ 



(59) 
(60) 
(61) 



The quantity ip(t p ) is defined in Section IIII Bl The 
asymptotic limits of Eqs. ([55]) through (|6"Tj) are given by 



Go- (a) _ n j 




T^a, a<10~ 2 , 



l 7r 3/2 -l/2 ) a >l 2, 
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a<0.1, 



U o < 10 _ a (62) 



For a values outside the ranges listed above the values of 
H and G were calculated numerically. Asymptotic ap- 
proximations for S a were used for t p < 10~ 2 and t p > 1. 
For t p values outside of these intervals the expression for 
S a was computed numerically. The asymptotic approx- 
imations for 5S K are valid, and have been used, for t p 
values in the intervals t p < 10~ 2 and t p > 10. 



D. Specific heat 

The calculation of the specific heat is much simpler 
than the calculation of the neutrino emissivities and of 
the thermal conductivity, since the specific heat does not 
involve scattering processes. In terms of the phonon sum 
used in Sections IIII Bl and IIII C[ the specific heat density 
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can be written as 

CV.solid = Swifc 



4*2 



, sinh' 



(2**) 



(63) 



The asymptotic low and high temperature expressions of 
Eq. (|63|) are given by 



3k B n b ( 1 



Cy.solid 



l 



t P > io 2 , 



12k B n b {^ + ^)tl, t P < 10- 



(64) 



For temperatures greater than the melting temperature 
we adopt Eq. (24) of [H| for the specific heat density, 



k B n b . 



CVAii 



l kB {} 1 io g r„ 



r < i 

i< r < r m , 



(65) 



where we use r m = 172 (as in Sections IIII Bl and IIII C| 
instead of the outdated value T m = 150 |44j |. 



IV. RESULTS AND DISCUSSION 



A. Contribution to the specific heat 
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FIG. 10. (Color online) Specific heat, cy, of mixed quark- 
hadron phase as a function of quark volume fraction, \i com- 
puted for different surface tensions, 7/, and temperatures, T. 
The underlying equation of state is HV1. The panel on the 
left hand side shows the contributions of the rare-phase blobs 
to the specific heat. The curves in the panel on the right 
hand side show the specific heat computed for a standard (no 
geometrical structures) quark-hadron gas. 



In this paper, we have calculated the specific heat 
stemming from rare phase blobs immersed in hadronic 
matter for four different parameter sets (see Table IIII[) 
and three different values (77 = 0.5, 1, 2) for the surface 
tension of rare phase blobs. It is intriguing to compare 
these results with the heat capacities of the hadronic and 
quark matter phases, weighted by their volume fractions. 
For this purpose we compute the specific heat of a Fermi 
gas of leptons and baryons from 

c v = % T \/ m ? + k h k ^ > ( 66 ) 

C V = ll T E m B k F,B , (67) 
B 

where m* B is effective in-medium mass of baryons (see 
Section QTJ). The specific heat of a quark gas is given by 



c q v = 0.6 x 10 20 ( y J T 9 ergs cnT 3 KT 1 . (68) 

The total specific heat of a mixture of quarks and hadrons 
follows from cy = (1 — x) c v + X c v where Cy = Cy + c l v 
and Cy = c v + c v . Figure [10] compares the different 
contributions to the specific heat with one another. As 
can be seen, the contribution of the rare-phase blobs to 
the specific heat (colored lines in Fig. [TU)) is typically 



several orders of magnitude smaller than the specific heat 
of a standard (no geometrical structures) quark-hadron 
gas (solid black lines in Fig. [TU]) . The small jumps at 
X = 0.5 in Figure [10] are due to the discontinuity of the 
rare phase blob mass, mi, at \ — 0-5 (see Fig. [?])• The 
parameter sets HV2, G300I and G300II lead to results 
very similar to those shown in Fig. [TU] and are therefore 
not shown separately. 



B. Neutrino Bremsstrahlung emissivity 

The neutrino Bremsstrahlung emissivities emerging 
from electron-phonon scattering and electron-lattice 
(Bragg diffraction) scattering have been computed from 
Eqs. (|33p for the four parameter sets of Table ITO]) . The 
surface tension of rare phase blobs in Eq. ([M]) has been 
varied again from r\ = 0.5, 1, to 2. Figures [TT1 and [T2l 
show the contributions of the rare phase blobs to the 
neutrino Bremsstrahlung emissivity as a function of the 
quark volume fraction. A range of representative temper- 
atures, from T = 10 7 K to 10 11 K has been chosen. For 
comparison, we show the contributions to the neutrino 
emissivity which comes from the modified Urea process in 
hadronic matter, c^h.mUj and in quark matter, c^q.mu 
[45l ] . Finally, we also show in these figures the emissivi- 
ties which correspond to nucleon Bremsstrahlung, c^h^r 
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FIG. 11. (Color online) Neutrino emissivities, e„, as a 
function of quark volume fraction, \, at different temper- 
atures, computed for parameter sets HV1 and HV2. The 
colored lines show the total contributions from electron-blob 
Bremsstrahlung for different surface tensions: r\ = 0.5 (red), 
r\ = 1 (green), r\ = 2 (blue). The black lines show the 
total contributions from modified nucleon and quark Urea 
processes (dashed line) and nucleon-nucleon and quark-quark 
Bremsstrahlung processes (dotted line). 

and quark Bremsstrahlung, e^c^BR @, H^. Their total 
contribution in quark-hybrid star matter is obtained, for 
a given quark volume fraction \, by f-v — X e Q + (1 — x) £ h- 
As can be seen from Figs. [TT1 and [T2l the neutrino emis- 
sivity from electron-blob Bremsstrahlung become com- 
parable to the emissivities of the modified Urea process 
(and other Bremsstrahlung processes) for temperatures 
T < 10 8 K. 

The Bremsstrahlung emissivities oscillate rapidly with 
X for T < 10 9 and \ > 0.5. This is due to Bragg diffrac- 
tion, given by the sum in the expression for L s \ (see Eq. 
(f5Tj) ). As mentioned in Section fill B[ the sum of L s \ con- 
sists only of a few terms if \ ^ 0.5. The oscillations are 
essentially due to the oscillating values of the individ- 
ual terms in Eq. (|5"Tj) . The oscillations are smoothed out 
when the number of terms is large. 

It also follows from Figs. fTTI and IT21 that the neutrino 
emissivity of rare phase blob Bremsstrahlung, at low tem- 
peratures but greater quark volume fractions (\ > 0.5), 
becomes sensitive to the choice of r\. This feature, how- 
ever, does not lead to large uncertainties in the total neu- 
trino emissivity, since the neutrino emissivities for this 
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FIG. 12. (Color online) Same as Fig. 1111 but for parame- 
ter sets G300I and G300II. The colored lines show the total 
contributions from electron-blob Bremsstrahlung for different 
surface tensions, r\ — 0.5 (red), r\ = 1 (green), 77 = 2 (blue). 
The black lines shown the total contributions from modified 
nucleon and quark Urea processes (dashed line) and nucleon- 
nucleon and quark-quark Bremsstrahlung processes (dotted 
line). 

temperature-density regime are dominated by the mod- 
ified Urea process and nucleon-nucleon and quark-quark 
Bremsstrahlung processes. 

C. Thermal conductivity 

Figures [T3l through [Tol show the thermal conductivities 
due to presence of rare phase blobs. For comparison, 
thermal conductivities in standard quark-hadron gas are 
shown too where the thermal conductivity for hadronic 
phase kh is adopted from [4fj| and that of quark phase 
kq is adopted from (47j | . The total thermal conductivity 
in standard quark-hadron gas without blobs is given by 

«=(—+/—) . (69) 

Besides scattering between electron and rare phase blobs 
calculated in i jlll C[ the geometric pattern will also alter 
the total thermal conductivity. The two contributions 
kh and kq from standard quark-hadron gas can be com- 
bined using an expression for the total effective thermal 
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FIG. 13. (Color online) Thermal conductivity, func- 
tion of quark volume fraction, x> for matter at a temperature 
of 10 7 K and 10 s K. The colored lines show the contribu- 
tions of electron-blobs scattering in the mixed quark-hadron 
phase to k, computed for parameter sets HV1 (solid) and HV2 
(dashed), for different surface tensions (see Eq. (|24[l 'l: r\ = 0.5 
(red), Tj = 1 (green), rj = 2 (blue). The thick black lines show 
thermal conductivities of a standard (no blobs) quark-hadron 
gas computed for HV1 (solid) and HV2 (dashed). The thin 
black lines show the effective bulk thermal conductivities de- 
fined in Eq. ((70)) in the presence of rare phase blobs embeded 
in dominant phase without contribution from electron-blob 
scattering for HV1 (dotted line) and HV2 (dash-dotted line). 

conductivity of spheres immersed in continuous matter 
of a different thermal conductivity (48[. In our case the 
bulk thermal conductivity of two phases can be written 

as 

Keff = Kl ( 1 - — ■ — J, (70) 

where m and K2 are the thermal conductivities of the 
dominant and the rare phase, respectively. 

The small jumps of thermal conductivities from 
electron-blob scattering (color lines in Figs. [13] through 
[in]) are due to discontinuities of rare phase blob mass m& 
(See Fig. [7]). The jumps of thermal conductivities con- 
tributed by embedding of rare phase blobs (thin black 
lines in Figs. [T3l through I16[) are due to unequal thermal 
conductivities of the two phases «i ^ K2 at x — 0.5. 

As can be seen from Figs. [13] through [T6l the total 
thermal conductivity k — (k~^ + K bl Q t>) -1 1S dominated 
by electron-blob scattering at T < 10 9 -?C This is partic- 
ularly the case for quark volume fractions close to \ = 
or x — 1 where the electron thermal conductivity from 
blob scattering can be as much as three (for x < 0.5) 
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FIG. 14. (Color online) Same as Fig. 1131 but for a temperature 
of T = 10 9 K and 10 10 K. 

to six orders of magnitude (x > 0.5) smaller than the 
contribution from the mixed quark-hadron phase. Physi- 
cally, this causes a blocking of the thermal flow through a 
mixed quark-hadron phase region, which could manifest 
itself in the thermal evolution of quark- hybrid stars. 

As mentioned in Section [II C[ the blob volume density 
rib — > as X ~ * or x ~ ^ 1- Since the Coulomb logarithm 
A K is finite for these limits but the coefficient oc Z e g is 
divergent, the thermal conductivity stemming from rare 
phase blob scattering will diverge on both ends of the 
quark-hadron phase, leading to a vanishing total thermal 
conductivities there. Since this occurs only very near 
the edges of the quark-hadron boundary (x < 10~ 2 ) this 
feature can not be seen in Fig. 1131 The small jumps in k 
near the edges of the quark-hadron phase for T = 10 9 K 
and 10 10 K (see Figs. ITU and HI)]) are due to melting. 

V. SUMMARY AND CONCLUSIONS 

Because of the competition between the Coulomb 
and the surface energies associated with the positively 
charged regions of nuclear matter and negatively charged 
regions of quark matter, the mixed phase may develop 
geometrical structures (e.g., blobs, rods, slabs), similarly 
to what is expected of the sub-nuclear liquid-gas phase 
transition. In this paper we explore the consequences of a 
Coulomb lattice made of rare phase blobs for the thermal 
and transport properties of neutron stars. The total spe- 
cific heat, cy, thermal conductivity, k, and electron-blob 
Bremsstrahlung neutrino emissivities, c^.br? are calcu- 
lated and compared with those of standard neutron star 
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FIG. 15. (Color online) Thermal conductivity, func- 
tion of quark volume fraction, x> for matter at a temperature 
of 10 7 K and 10 s K. The colored lines show the contribu- 
tions of electron-blobs scattering in the mixed quark-hadron 
phase to k, computed for parameter sets HV1 (solid) and HV2 
(dashed), for different surface tensions (see Eq. ([24} ) : r\ = 0.5 
(red), Tj = 1 (green), 77 = 2 (blue). The thick black lines show 
thermal conductivities of a standard (no blobs) quark-hadron 
gas computed for G300I (solid) and G300II (dashed). The 
thin black lines show the effective bulk thermal conductivities 
defined in Eq. I|70|l in the presence of rare phase blobs embe- 
ded in dominant phase without contribution from electron- 
blob scattering for G300I (solid) and G300II (dashed). 

matter. To carry out this project, we have adopted, and 
expanded on, methods of earlier works on the transport 
properties of neutron stars [7, 23]. The sizes of, and spac- 
ings between, rare phase blobs are calculated using the 
Wigner-Seitz approximation Q. The equations of state 
used in this study are computed for a standard non-linear 
nuclear Lagrangian, and the associated equations of mo- 
tion for the baryon and meson fields are solved in the 
relativistic mean-field approximation. Quark matter has 
been modeled in the framework of the MIT bag model. 
Four different parameter sets (HV1, HV2, G300I, G300II) 
have been used to model the composition of neutron star 
matter containing a mixed phase of quarks and hadrons 
(quark- hybrid matter). 

The results discussed in Section llVl show that the con- 
tribution of rare phase blobs in the mixed phase to the 
specific heat is negligible compared to the specific heat of 
a quark-hadron gas. This is very different for the trans- 
port properties. For low temperature T < 10 s K the neu- 
trino emissivity from electron-blob Bremsstrahlung scat- 
tering is at least as important as the total contribution 
from other Bremsstrahlung processes (such as nucleon- 



FIG. 16. (Color online) Same as Fig. 1151 but for a temperature 
of T = 10 9 K and 10 10 K. 



nucleon and quark-quark Bremsstrahlung) and modified 
nucleon and quark Urea processes (see Figs. ITTI and IT2"]) . 
It is also worth noting that the scattering of degener- 
ate electrons off rare phase blobs in the mixed phase re- 
gion lowers the thermal conductivity by several orders of 
magnitude compared to a quark-hadron phase without 
geometric patterns (see Figs . [131 through [T6|) . This may 
lead to significant changes in the thermal evolution of the 
neutron stars containing solid quark-hadron cores, which 
will be part of a future study. 

In summary, our study has shown that the presence 
of rare phase blobs in dense neutron star matter may 
have very important consequences for the total neutrino 
emissivity and thermal conductivity of such matter. The 
implications of this for the thermal evolution of neutron 
stars need to be explored in future studies. To accom- 
plish this we intend on performing two dimensional cool- 
ing simulations, in which rotation and a dynamic compo- 
sition might be accounted for [49|, [5(| ■ In this connection 
we refer to the recent study of Noda et al. [5l| , who sug- 
gested that the rapid cooling of the neutron star in Cas- 
siopeia A can be explained by the existence of a mixed 
quark-hadron phase in the center of this object. 
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